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Abstract 

In this article, we briefly review recent progress related to extracting the quark-gluon plasma 
(QGP) specific shear viscosity from the flow data measured at Relativistic Heavy-Ion Collider 
(RHIC) and the Large Hadron Collider (LHC). 



1. Introduction 

Heavy ion experiments at the Relativistic Heavy- Ion Collider (RHIC) at BNL and the Large 
Hadron Collider (LHC) at CERN have accumulated strong evidences for the creation of the quark 
gluon plasma (QGP) that is strongly coupled and behaves like an almost perfect liquid. The 
transport properties of the QGP fluid is a currently hot topic. Due to the difficulties of the first 
principle calculations, it is desirable to extract them from experimental data. In this proceeding, 
we will review recent progress related to extracting the QGP shear viscosity from the flow data 
at RHIC and the LHC. 



2. Hydrodynamics and hybrid models 

Viscous hydrodynamics is an useful tool to study the viscous effects on the QGP fireball 
evolution and final observables for relativistic heavy-ion collisions at RHIC and the LHC ener- 
gies. During the past years, several groups have independently developed (2+l)-d [Q], 0. 0] and 
(3+l)-d |4|,|5|] viscous hydrodynamic codes with/without longitudinal boost invariance. A recent 
comparison between 2+1-d and 3+1-d viscous hydrodynamics shows that the realistic longitudi- 
nal expansion from the 3+1-d code only slightly affects the flow at mid-rapidity Jfl]. Therefore, 
one can safely implement the 2+1-d code to investigate the soft physics at mid-rapidity. 

With the efforts from different viscous hydrodynamic groups, flow has been widely accepted 
as the key observable to extract the QGP viscosity. However, the hadronic chemical composi- 
tion and non-equilibrium kinetics also play important roles for flow development in the hadronic 
phase, bringing additional uncertainties for the extracted QGP viscosity. For a better descrip- 
tion of the hadronic matter, one further developed the viscous hydrodynamics + hadron cascade 
hybrid approach that connects the hydrodynamic evolution for a viscous QGP fluid to the micro- 
scopic description for the re-scatterting and decoupling of the late hadronic matter. Such hybrid 
approaches include VISHNU hybrid model [7] developed in 2010 based on 2+1-d viscous hydro- 
dynamics and UrQMD hadron cascade (see Sec 4 for a summary of the main results), the McGill 
code developed in this year based on 3 + 1 -d viscous hydrodynamics and UrQMD (the new results 
can be found in Ref. ]8|]) and the Livermore hybrid code (see Sec. 7 and references therein). 
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Single shot vs. event-by-event simulations: Many hydrodynamic calculations concentrate on 
single-shot simulations for computing efficiency. Recently, several groups extended single-shot 
hydrodynamic simulations to event-by-event ones, making it possible to investigate the hydrody- 
namic response of the initial state fluctuations J3,l2, 13. Single shot simulations with smooth 
initial conditions obtained by averaging over large numbers of fluctuating events after rotation 
to align the "event planes" are, in principle, a computationally efficient short-cut to investigate 
initial-state fluctuation effects on even and odd flow coefficients. However, a detailed compar- 
ison between event-by-event and single-shot hydrodynamics shows 0(10%) deviations for the 
elliptic and triangular flows, especially when the QGP viscosity is small ll ill . Furthermore, 
the higher order flow harmonics V4, V5 and \% can not be reliably computed within a single-shot 
hydrodynamic approach due to the mode-coupling effects. With the appearance of new flow 
measurements (event-by-event v„ distribution, event plan correlation between flow angles, etc., 
please refer to Sec. 7 for details) and the ever increasing computing and storage capabilities, the 
transition to event-by-event simulations are inevitable. 



3. The QGP shear viscosity from the elliptic flow data - an early attempt 

The elliptic flow and higher order flow harmonics are important observables for the bulk 
matter created in relativistic heavy ion collisions. It was found that the elliptic flow V2 and 
triangular flow V3 are very sensitive to the shear viscosity. Even the conjectured lowest value 
rj/s — l/(47r) from the AdS/CFT correspondence leads to a significant suppression of V2 and 
v 3 01 SB 01- The bulk viscosity also suppress the development of flow anisotropy. Its effects 
are much smaller than the shear viscous ones due to the critical slowing down near the phase 
transition. Therefore, one can extract the QGP shear viscosity from flow data without large 



contaminations from the bulk viscosity [ 12] 



The first attempt to extract the QGP viscosity from the elliptic flow data, using 2+1-d viscous 
hydrodynamics, was done by Luzum and Romatschke around 2008 J2]. They implemented two 
initial conditions from optical Glauber and KLN models and found that the ~30% uncertainties 
in the initial eccentricity lead to ~30% uncertainties for the calulated elliptic flow, which trans- 
lates into ~100% uncertainties for the extracted value of the QGP shear viscosity. The effects 
neglected in this work are the off-equilibrium kinetics (or so-called highly viscous hadronic ef- 
fects) and the partially chemically equilibrated nature of the hadronic phase, and the initial state 
fluctuations. Around that time, the effects from bulk viscosity was unknown. After making gen- 
erous estimations for all of these uncertainties, it appears that the averaged specific QGP shear 
viscosity, cannot exceed the following conservative upper limit 131: 



< 5 x — . 
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4. The QGP shear viscosity at RHIC and LHC energies - the current status 

4.1. The QGP viscosity at RHIC energies 

With VISHNU hybrid model becoming available and the effects from hadronic matter being in 
control, we were ready to extract the QGP shear viscosity from the flow data with reliable uncer- 
tainty estimates. Ref. 01411 proposes to extract the QGP viscosity from the integrated V2 data for all 
charged hadrons since it is most directly related to the fluid momentum anisotropy and less sen- 
sitive to the details of VISHNU calculations, such as the form of the non-equilibrium distribution 
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Figure 1 : (Color online) eccentricity-scaled elliptic flow as a function of final multiplicity per area. 



function 5f, the bulk viscosity, initial flow, etc QUI. it was found that the eccentricity-scaled 
elliptic flow Vj'Vs.v as a function of charged multiplicity per unit overlap area (l/S)dN c h/dy is 
universal that depends only on the QGP shear viscosity but not on the initialization models 11411 . 
It thus desirable to use these curves to extract the QGP viscosity. Fig. 1 shows a comparison 
of the theoretical and experimental v c J' /s x - (l/S)dN c h/dy curves, where the left and right pan- 
els correspond to MC-Glauber and MC-KLN initialization models. The theoretical lines are 
calculated from single shot VISHNU simulations with different (tj/s)qcp as inputs, using event- 
averaged smooth initial conditions generated from these two models. The experimental data are 
the corrected elliptic flow with non-flow and fluctuation effects removed ifTetl . Due to the -20% 
larger ellipticity of the MC-KLN fireballs, the magnitude of \" p /s x differs between the two mod- 
els. As a result, the extracted value of (j]/s)qcp from these two panels changes by a factor of 
2. Taken the main uncertainties from these two initial conditions, we found that the averaged 



specific shear viscosity of the QGP created at top RHIC energies is 11411 : 



1<2 
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Remarks: Small bulk viscous effects and proper event-by-event hydrodynamical evolution of 
fluctuating initial conditions slightly reduce the integrated V2. while early flow slightly increase 
it. Although they should be studied in more quantitative detail, we expect the total uncertainty 
band translated to the above extracted value of QGP shear viscosity may slightly shift after can- 



celations. For more details on these uncertainty estimates, please refer to 1171. 

With (t]/s) qgp - (l/An) for MC-Glauber and (i]/s) QGP {2/ An) for MC-KLN extracted 
from the pj integrated vn for all charged hadrons, VISHNU yields an excellent description of 
the /?7-spectra and differential elliptic flow V2(pr) (for pj < 2GeV) for all charged hadrons 
and identified hadrons at different centrality bins measured in 200 A GeV Au+Au collisions at 



RHIC 11511 . showing the power of the VISHNU hybrid model and the robustness of the above 



extraction of the QGP shear viscosity. 



4.2. Recent results at the LHC 

Elliptic flow and the QGP shear viscosity at the LHC: After extrapolating our calculation to 
the 2.76 A TeV Pb+Pb collisions at the LHC, we found that the 0]/s)q GP extracted from RHIC 
slightly over-predicts the LHC elliptic flow data. After increasing (t]/s)q C p from ~ 2/ (An) to 
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~ 2.5 1 (An) (for MC-KLNfl VISHNU gives a better description of the integrated and differential 
elliptic flow data for all charged hadrons at the LHC IU9I1 . However, this dost not necessarily 
mean that the QGP turns more viscous at higher temperatures reached by LHC since vi becomes 
sensitive to the initial shear stress tensor at the LHC and one needs to further investigate the 
temperature-dependent QGP viscosity (please refer to Sec. 6 for detailed discussions). 

The ALICE collaborations also measured the differential elliptic flow v-l(Pt) for identified 
hadrons from central to peripheral collisions. Although purely viscous hydrodynamic simula- 
tions makes a nice description of the pion and kaon data at various centralities, it over-predicts 
the proton V2(pr) from most central to semi-central collisions |2^, 21]. This problem was fixed 



by the later VISHNU hybrid model simulations. With the microscopic rescattering of the hadronic 
matter, the generations of radial and elliptic flow for protons are re-balanced. Besides keeping 
the nice descriptions of the pion and kaon data, VISHNU gives an improved description of the 
proton elliptic flow data at various centrality bins, especially for the most central to semi-central 
collisions 12111 

A simultaneous fitting of the elliptic and triangular flow data: The QGP shear viscosity sup- 
presses the development of both elliptic flow v>2 and triangular flow V3 |H S S H Ell- While 
the initial eccentricities £2 differ by 0(20%) between MC-KLN and MC-Glauber models, the 
triangular deformation £3 are almost identical. As a result, V3 is much less sensitive to these two 



initializations when compared with V2 II Hi . Using pure viscous hydrodynamics, a simultaneous 
fitting of the elliptic and triangular flow data indicates that these recent LHC data favors a small 
value of QGP shear viscosity around 1 /(An). With that value, the centrality dependent integrated 
elliptic and triangular flow are well described by the MC-Glauber initial conditions, while MC- 
KLN initial conditions obviously under-predict the triangular flow data with the QGP viscosity 
tuned to fit the elliptic flow (22]. Although more comprehensive investigations from the comput- 
ing expensive VISHNU hybrid model are needed in the near future, the hadronic afterburner will 
not change the main conclusion since both elliptic and triangular flow almost reach saturation in 
the QGP phase at the LHC energies lEI . 

The proton puzzle at the LHC: Recently, ALICE collaborations measured the p/n~ ratio in 
2.76 A TeV Pb+Pb collisions at the LHC 12411 . which obviously disagrees with the value observed 
by STAR in 200 A GeV Au+Au collisions at RHIC, but agrees with the one from PHENIX. It 
was found that p/n~ ratio is sensitive to the chemical freeze-out temperature T c i tem . in theoretical 
calulations. Correspondingly, hydrodynamic simulations with T c ] mlL = 165 MeV (chosen the 
value measured at RHIC) over-predict the p spectra at the LHC [20]. Our past VISHNU simula- 
tions (that accidentally tuned off the baryon anti-baryon B - B annihilation channels) with the 
switching temperature T sw - T c hem. = 165 MeV also over-predict the p spectra at the LHC lfl9ll . 
To solve this problem, some groups proposed to use a lower chemical freeze-out temperature to 
fit the data. However, new VISHNU simulations with the B - B annihilation channels shows that 
such annihilations reduce the proton and anti-proton multiplicity by about 30% in central col- 
lisions and by about 15% in peripheral collisions, leading a dramatically improved description 
of the proton and anti -proton spectra measured at the LHC 12511 . With that B - B annihilation, 
the descriptions of the proton and anti-proton final-multiplicity and spectra at RHIC (measured 



by PHENIX) are also improved 12511 . Although the detailed work is still on-going, this helps to 



explain the theoretical over-prediction of the proton, anti-proton spectra, and p/n ratio at the 



For computing efficiency, we only investigate the case with MC-KLN initializations in details. The trend on how 
(*1/ s )qgp changes from RHIC to LHC energies is similar for MC-Glauber and MC-KLN initializations flUl . 

4 



LHC with the chemical freeze-out temperature set to 165 MeV. 



5. Recent developments on initialization models 

The MC-Glauber and MC-KLN models, commonly used by many groups, construct the ini- 
tial state fluctuations through the positions fluctuations of the nucleons inside the colliding nuclei. 
Recently, several groups further investigate the quantum fluctuations of color charges. Combing 
the classical Yang-Mill's approach for the Glasma field with the impact parameter dependent 
saturation model (called IP-Glama model), Schenke, Tribedy and Venugopalan studied the color 
charge fluctuations and found that such effects moderately change the initial eccentricities com- 



pared with the one from traditional MC-KLN model. Ref. 12711 investigated the transverse corre- 
lations for energy density fluctuations within the framework of Color Glass Condensate (CGC). 
Based on this work, the OSU group constructed a new Monte-Carlo initialization generator to 
produce initial energy density profiles with correlated fluctuations and found a small increase 
of 62 - £5 In Ref- i29ll . the fluctuations of the initial gluon productior0 are investigated 

accounting to a negative binomial distribution within the kj factorization approach of CGC. It 
gives an initial eccentricity £2 close to the one from traditional MC-KLN model, but obvious 
larger £3 - £5 than the MC-KLN one. The larger higher order eccentricity leads to larger value of 
theoretical higher order flow harmonics, which requires a significantly larger value of the QGP 
shear viscosity to fit the corresponding experimental data. This brings additional uncertainties 
for our extracted (77/ s)qqp from the elliptic flow data alone. A joint investigation of elliptic, tri- 
angular and higher order flow harmonics can help to further constrain the initialization models 
and yield a more accurate value for (77/ s)qqp. 

Hydrodynamic initial conditions can also be provided by dynamical models, such as URQMD 



113011 , EPOS 1311], AMPT Q32D, IP-Glama 126J], etc., which try to account for the pre-equilibrium 
dynamics before thermalization. Generally, pre-equilibrium dynamics reduces the initial eccen- 
tricities £2 and £3 and contributes radial and anisotropic flow. It was found that the additional 
fluctuations from the energy deposition in UrQMD slightly increase £2 and £3 ll30Tl and the early 
longitudinal flow fluctuations from AMPT slightly reduce V2 Except for IP-Glasma, most 
of these dynamical models are matched to ideal hydrodynamics for further evolution. How they 
quantitatively influence the extracted value of the QGP shear viscosity is still unknown and needs 
additional work. 



6. Temperature-dependent (ti/s)q GP (T) 

Besides further increasing the precision of the extracted "averaged" specific QGP shear vis- 
cosity, it is important to quantitatively determine the temperature-dependent (77/ s)qcp{T) to fur- 
ther explore the fundamental properties of the QCD matter at high temperatures. Past research 
showed that one needs a slightly larger constant 77/s to fit the LHC flow data than the one used 
for RHIC Jl9l 351. which demonstrates that it is feasible to extract a temperature-dependent 



(jiI s)qcp(T) from the new measured flow data at different collision energies. One of the crucial 
issues is to extract the initial temperature since it can not be directly measured in experiments. 
The new developed massive data evaluating technique [34] make it possible to simultaneous ex- 
tract the QGP shear viscosity, the initial temperature and other parameters in initial conditions 



2 For other models including multiplicity fluctuations for soft particle productions, please refer to Ref. ( 10, 30, 32, 33]. 
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from multiple data sets. The error bands of these extracted values are largely controlled by the 
precision of the measured experimental data. It is thus desirable to further reduce the systematic 
and statistical uncertainties of related observables (such as final multiplicities and the pj spectra, 
HBT radii, flow harmonics v„, event plan correlations of flow angles, v„ distributions, etc., please 
also refer to Sec. 7 for details) for a quantitatively determination of (rjl s)qcp(T). 

7. Other related current developments 

Event plane correlations for flow angles: Recently ATLAS collaborations measured the event 
plan correlations between the flow angles associated with higher order flow harmonics in 2.76 A 
TeV Pb+Pb collisions 13611 . Event-by-event hydrodynamic simulations show that the correlation 
strength is sensitive to both the initial conditions and the QGP shear viscosity l37ll . Although 
more detailed studies are needed for a quantitative description of the data, it is impressive that 
e-by-e hydrodynamics correctly reproduced the qualitative features of the centrality dependent 
correlations between different flow angles, providing strong support for the fluid dynamic de- 
scription of the fireball evolution. 

Event-by-event distribution of v n : ATLAS also measured the event-by-event distributions of 
vn - v>4 at varies centralities in 2.76 A TeV Pb+Pb collisions 13811 . The event-by-event hydrody- 
namic simulations with MC-Glauber and MC-KLN initializations show that none of these two 
initializations works for the distributions for all v„ l39ll . In contrast, IP-Glasma gives the s n 
distributions that largely overlap with these measured v„ distributions. After the hydrodynamic 
evolution, the descriptions of the data are further improved ll3~5ll . A future study of these flow dis- 
tributions with different initializations as inputs may help us to constrain the initialization models 
and to understand the sources of fluctuations. 

Higher order flow harmonics in ultra-central collisions: In Ref. l40ll . Luzum proposed to ex- 
tract the QGP viscosity from a simultaneous fitting of all measured v„ in the ultra-central colli- 
sions since the fluctuation effects are dominated and the geometry effects are suppressed there. 
However, viscous hydrodynamic simulations from the OSU group, with different rj/s as input, 
show that neither the MC-KLN nor the MC-Glauber model appears to be able to simultaneously 
describe all measured v„ from CMS I4l[|. It will be interesting to see a comparison of the other- 
wise very successful IP-Glasma model 13511 with these ultra-central Pb-Pb data in the near future. 

Systematic^ 2 fitting of the experimental data: Using a hybrid model that connect UVH2+1D 
viscous hydrodynamics with UrQMD hadron cascade model, Soltz and his collaborators made a 
systematic x 2 evaluation of the pion spectra, elliptic flow and HBT radii measured by STAR 
and PHENIX in 200 A GeV Au+Au collisions 13411 . The evaluation is performed in a two- 
dimensional parameter space for initial temperature and the QGP shear viscosity. For each in- 
putting initial conditions (with/without pre-equilibrium flow), the fitted value for them is give by 
the lowest value of the sum of x 1 for these three observables. Although this work only concen- 
trate on the data sets in the 0-20% centrality bin and has lots of room for future improvement, it 
demonstrates a systematic evaluation technique for constraining the multiple theoretical param- 
eters from the multiple experimental data sets for relativistic heavy ion collisions. 

8. Summary and concluding remarks 

In summary, our field has evolved quickly during the past few years. Using pure viscous 
hydrodynamics, an early extraction of the QGP shear viscosity from elliptic flow data suggested 
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Figure 2: (Color online) Time line for the extracted QGP specific shear viscosity with uncertainty estimations. Grey 
bands are for RHIC (200 A GeV Au+Au collisions) and dark green band are for LHC (2.76 A TeV Pb+Pb collisions). 



{iiI s )qgp < 5/(4tt) I0L with the main uncertainties coming from the undetermined initial con- 
ditions, an improper description of hadronic chemical composition and non-equilibrium kinetics 
of the late hadronic stage. With the hadronic afterburner that realistically describes these two 
hadronic effects, the VISHNU hybrid model fits the centrality-dependent integrated elliptic flow 
data with MC-Glauber and MC-KLN initial conditions and gives 1 /(An) < (t]/s)q GP < 2.5 1 {An), 
where the uncertainty is now entirely dominated by the different initial eccentricities from these 
two models [14]. One also found that the averaged QGP shear viscosity at the LHC is close to (or 
slightly larger than) the one at RHIC lfl9ll . Using pure viscous hydrodynamics, a simultaneous 
fit of the elliptic and triangular flow at 2.76 A TeV Pb+Pb collisions strongly indicated that these 
data prefer an even smaller value of the specific QGP viscosity around l/(47r). With that value, 
the MC-Glauber initial conditions can fit the centrality dependent elliptic and triangular flow 
data, whereas MC-KLN initial conditions fail to do so 12211 . Since both elliptic and triangular 
flow almost reach saturation in the QGP phase at the LHC, future sophisticated VISHNU calcula- 
tions with the hadron afterburner are not expected to change these conclusions. The MC-Glauber 
and MC-KLN initial condition models used by many hydrodynamic simulations (including our 
own) construct the initial state fluctuations from nucleon position fluctuations inside the collid- 
ing nuclei. Recently, several groups further studied other sources of quantum fluctuations in the 
initial state, such as color charge fluctuations, multiplicity fluctuations, initial flow fluctuations, 
etc. [26, 27, 28, 2§l 32]. Initially this led to large additional uncertainties for the extracted QGP 
shear viscosity. A recent fit of v„ and event-by-event v„ distributions with 3+1-d viscous hy- 
drodynamics using IP-Glasma initial conditions 0511 dramatically reduced these uncertainties, 
by eliminating both the MC-Glauber and MC-KLN models as viable models for the initial-state 
fluctuations 051, 13911 . The combined analysis of all measured v„ (« = 1, ... ,5) gives tj/s ~ 0.20 
for LHC collisions and the smaller value tj/s ~ 0.12 for RHIC collisions 0511 . Although detailed 
studies of initialization models are still on-going, additional systematic anisotropic flow stud- 
ies (including the centrality dependence of all integrated and differential flow harmonics for all 
charged and identified hadrons, event plane correlations between flow angles, flow distributions, 
etc.) are expected to strongly constrain the initial conditions and to further narrow the error band 
for the extracted QGP viscosity. Progress is happening quickly, and exciting new results are 
expected to be available by the next Quark Matter meeting. 
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